The effects of spatial constraints on the evolution of weighted complex networks 
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Motivated by the empirical analysis of the air transportation system, we define a network model 
that includes geographical attributes along with topological and weight (traffic) properties. The 
introduction of geographical attributes is made by constraining the network in real space. Inter- 
estingly, the inclusion of geometrical features induces non-trivial correlations between the weights, 
the connectivity pattern and the actual spatial distances of vertices. The model also recovers the 
emergence of anomalous fluctuations in the betweenness-degree correlation function as first observed 
by Guimera and Amaral [Eur. Phys. J. B 38, 381 (2004)]. The presented results suggest that the 
interplay between weight dynamics and spatial constraints is a key ingredient in order to understand 
the formation of real-world weighted networks. 
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I. INTRODUCTION 



The empirical evidence coming from studies on sys- 
tems belonging to areas as diverse as social sciences, bi- 
ology and computer science have shown that the usual 
paradigm of random graphs is often not well suited to 
describe real world networks 0,0, SB- In particular, in 
a wide range of networks the occurrence of vertices with 
a very large degree (number of links to other vertices) is 
very likely. The presence of these "hubs" often goes along 
with very large degree fluctuations. The large topological 
heterogeneity associated to these features is statistically 
expressed by the presence of heavy-tailed degree distri- 
butions with diverging variance that have a very strong 
impact on the networks' physical properties such as re- 
silience and vulnerability, or the propagation of pathogen 
agents HI 0,0. 

The purely topological definition of networks, however, 
misses important attributes which are frequently encoun- 
tered in real- world networks. In the first instance, net- 
works are far from boolean structure and are better rep- 
resented as weighted graphs with the intensity of links 
that may vary over many orders of magnitude. Indeed, 
in many graphs ranging from food-webs to metabolic net- 
works, large variations of the link intensities are empir- 
ically observed H E E E E H EI Notably, the 
statistical properties of weights indicate non-trivial cor- 
relations and association with topological quantities EH • 
Finally, the correlation between weights on different links 
is at the origin of the existence of pathways which are 
particularly important in metabolic networks for exam- 
ple |l5j . Another important element of many real net- 
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works is their embedding in the real space. For instance, 
most people have their friends and relatives in their 
neighborhood, transportation networks depend on dis- 
tance, and many communication networks devices have 
short radio range El E E E HI A particularly 
important example of such a "spatial" network is the In- 
ternet which is a set or routers linked by physical cables 
with different lengths and latency times [j, |21| . An anal- 
ogous situation is faced in the air transportation network 
with routes covering very different distances. The length 
of the link is a very important quantity usually associated 
with an intrinsic cost in the establishment of the connec- 
tion. If the cost of a long-range link is high, most of 
the connections starting from a given node will go to the 
closest neighbors in the embedding space. Long-range 
links, on the other hand, correspond usually to connec- 
tions towards already well-connected nodes (hubs). This 
seems natural in the case of the air transportation net- 
work for instance: short connections go to small airports 
while long distance flights are directed preferentially to- 
wards large airports (i.e. well connected nodes). It is 
therefore natural to find that spatial constraints can have 
important consequences on the topology of the resulting 
network 

Recently, the raising interest on the dynamics and 
function of complex networks has fostered studies go- 
ing beyond the simple topological structure. In partic- 
ular, models of complex networks in which the diversity 
of weights is taken into account have been formulated 
having in mind growing networks where the dynamics is 
driven by the intensity of the weights along with a re- 
inforcement mechanism [3 HH E E3 ■ Other models 
have focused on more geometrical mechanisms or some- 
how different dynamical rules H E, E EH- These 
models, however, are not able to reproduce all of the fea- 
tures observed in real world networks. For instance, the 
anomalous centrality fluctuations observed in Ref. E| do 
not find a rationalization in models based only on topol- 
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ogy and weight properties. On the other hand, some of 
these interesting and non-trivial features can result from 
the introduction of spatial attributes in the models' con- 
struction [2^. In this article, we discuss the interplay 
of the three aforementioned ingredients (heterogeneous 
topology, weights and spatial constraints) in a model 
of growing network combining these ingredients at once. 
The proposed model is obtained as the embedding of the 
weighted growing network introduced in |24| in a two- 
dimensional geometrical space. Spatial constraints are 
translated into a preference for short links, and combined 
with the coupling between the evolution of the network 
and the dynamical rearrangement of the weights. This 
mechanism naturally leads to the appearance of many 
features observed in real- world networks, in particular 
the non-linear correlations between weights and topology, 
and the large fluctuations of the betweenness centrality. 

The paper is organized as follows. In section [HI we 
briefly review some important empirical results of the 
North- American airline network, highlighting the most 
salient effects on space on different quantities. Sections 
IIIII and IIVI are devoted to the presentation and to the 
study of the spatial weighted model, stressing the effect 
of the spatial embedding and constraints on the proper- 
ties of the resulting network. In section [V] we present 
a summary of the results and conclusions about large 
network modeling. 



II. A CASE STUDY: SPACE, TOPOLOGY AND 
TRAFFIC IN THE NORTH AMERICAN 
AIRLINE NETWORK 

The airline transportation infrastructure is a 
paramount example of large scale network which 
can be represented as a complex weighted graph: the 
airports are the vertices of the graph and the links 
represent the presence of direct flight connections 
among them. The weight on each link is the number of 
maximum passengers on the corresponding connection. 
The characteristics of the world-wide air-transportation 
network using the International Air Transportation 
Association (IATA) database [34j have been presented 
in 0| . The network is made of N = 3880 vertices 
and E = 18810 edges and shows both small-world 
and scale-free properties as also confirmed in different 
datasets and analyses 0, 0,|23|,|32. In particular, the 
average shortest path length, measured as the average 
number of edges separating any two nodes in the network 
shows the value (£) = 4.37, very small compared to the 
network size N. The degree distribution, takes the form 
P(k) = k~ 7 f(k/k x ), where 7 ~ 2.0 and f(k/k x ) is an 
exponential cut-off function. The degree distribution is 
therefore heavy-tailed with a cut-off that finds its origin 
in the physical constraints on the maximum number of 
connections that a single airport can handle |23l . 133 . l33| . 
The airport connection graph is therefore a clear example 
of small-world network showing a heavy-tailed degree 




FIG. 1: Cumulative degree distribution P c (k) for the North- 
American network. The straight line indicates a power-law 
decay with exponent 7 — 1 = 0.9. 

distribution and heterogeneous topological properties. 

The world-wide airline network necessarily mixes dif- 
ferent effects. In particular there are clearly two different 
scales, global (intercontinental) and domestic. The inter- 
continental scale defines two different groups of travel 
distances and for the statistical consistency we elimi- 
nate this specific geographical constraint by focusing on 
a single continental case. Namely, in the following we 
will consider the North- American network constituted of 
N = 935 vertices with an average degree (k) w 8.4 and an 
average shortest path <Rii The statistical topological 
properties of the North American network are consistent 
with the world-wide one as we will see in the forthcoming 
analysis. 

A. Topology and weights 

The North American network presents a degree distri- 
bution statistically consistent with the world-wide airline 
network. Indeed, we also observe (Fig. in this case a 
power-law behavior on almost two orders of magnitude 
followed by a cut-off indicating the maximum number of 
connections possible due to limited airport capacity and 
to the size of the network considered. 

The existence of a broad degree distribution signals a 
strong heterogeneity of the network at the topological 
level which exists also at the weight level. A first indica- 
tion on the weight heterogeneity is given by the s tudy of 
the weight strength of a node i as defined by [Tj, |35| 

sf = J2 *>u ■ (i) 

jev(i) 

where the sum runs over the set V(i) of neighbors of 
i. The strength generalizes the degree to weighted net- 
works and in the case of the air transportation network 
quantifies the traffic of passengers handled by any given 
airport. This quantity obviously depends on the degree 
k and increases (linearly) with k in the case of random 
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FIG. 2: Weight and distance strengths versus degree for the 
North- American network. The dashed lines correspond to the 
power- laws /3d ~ 1.4 and /3 m ~ 1.7. 



FIG. 3: Distribution of distances (in kms) between airports 
linked by a direct connection for the North- American network. 
The straight line indicates an exponential decay with scale of 
order 1000 km. 



uncorrelated weights of average (w) . A relation between 
the average strength s w (k) of nodes of degree k of the 
form 

s w = Ak p ™ , (2) 

with an exponent j3 w > 1, or [3 W = 1 but A ^ (w) is 
then the signature of non-trivial statistical correlations 
between weights and topology. This is indeed what we 
observe in the North- American air transportation net- 
work with f3 w ~ 1.7 (Fig. EJl. 

B. Spatial analysis 



with j3 d ~ 1.4 (Fig. |2J). This result shows the presence 
of important correlations between topology and geogra- 
phy. Indeed, the fact that the exponents appearing in 
the relations J5J and (@J are larger than one have differ- 
ent meanings. While Eq. J5J means that larger airports 
have connections with larger traffic, implies that they 
have also farther-reaching connections. In other terms, 
the traffic (and the distance) per connection is not con- 
stant but increases with k. As intuitively expected, the 
airline network is an example of a very heterogeneous 
network where the hubs have at the same time large con- 
nectivities, large weight (traffic) and long-distance con- 
nections [l4|, related by super-linear scaling relations. 



The spatial attributes of the North American airport 
network are embodied in the physical spatial distance, 
measured in kilometers or miles, characterizing each con- 
nection. Fig. |3 displays the histogram of the distances 
of the direct flights. These distances correspond to Eu- 
clidean measures of the links and clearly show a fast de- 
caying behavior reasonably fitted by an exponential. The 
exponential fit gives a value for a typical scale of the or- 
der 1000 kms. The origin of the finite scale can be traced 
back to the existence of physical and economical restric- 
tions on airline planning in a continental setting. 

Since space is an important parameter in this network, 
another interesting quantity is the distance strength of i 

S f= £ dy (3) 

j'ev(i) 

where dy is the Euclidean distance. This quantity gives 
the cumulated distances of all the connections from (or 
to) the considered airport. Similarly to the usual weight 
strength, uncorrelated random connections would lead 
to a linear behavior of s d (k) oc k while we observe in the 
North- American network a power law behavior 

s d (k) ~ kfi* (4) 



C. Assortativity and Clustering 

A complete characterization of the network structure 
must take into account the level of degree correlations 
and clustering present in the network. Correlations can 
be probed by inspecting the average degree of the nearest 
neighbors of a vertex i 

knn,i = ^ ^ kj . (5) 

1 jev(i) 

Averaging this quantity over nodes with same degree k 
leads to a convenient measure to inve stig ate the behavior 
of the degree correlation function 0, l36| | 

knn{k*) — ~J^~ ^ ^ knn,i , (6) 
i/ki—k 

where iVj, is the number of nodes of degree k. This quan- 
tity 10 is related to the correlations between the degrees 
of connected vertices since on average it can be expressed 
as 

fc nn (fc) = £>'P(fc'|fc) . (7) 

k> 
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where P(k'\k) is the conditional probability that a given 
vertex with degree k is linked to a vertex of degree k' . 
If the degrees of neighboring vertices are uncorrelated, 
P(k'\k) is only a function of k' and thus k nn (k) is a con- 
stant. When correlations are present, two main classes 
of possible correlations have been identified: Assortative 
behavior if k nn (k) increases with k, which indicates that 
large degree vertices are preferentially connected with 
other large degree vertices; and dis assortative if k nn (k) 
decreases with k |37j . The weighted generalization of the 
above quantity, the affinity, reads as 0] 



k w = — V 

nnA / . 

Si '-^ 

iev(i) 



i^ij kj . 



(8) 



In this case, we perform a local weighted average of 
the nearest neighbor degree according to the normalized 
weight of the connecting edges, Wij/si. This definition 
implies that k™ n { > k nn ^ if the edges with the larger 
weights are pointing to the neighbors with larger degree 
and k™ n i < k nn ,i m the opposite case. The k™ n i thus 
measures the effective affinity to connect with high or 
low degree neighbors according to the magnitude of the 
actual interactions. As well, the behavior of k™ n (k), i.e. 
the affinity of vertices of degree k, marks the weighted 
assortative or disassortative properties considering the 
actual interactions among the system's elements. 

Information on the local connectedness is provided by 
the clustering coefficient Cj defined for any vertex i as 
the fraction of connected neighbors of i [33]. The aver- 
age clustering coefficient C — N~ 1 J2 i a thus expresses 
the statistical level of cohesiveness measuring the global 
density of interconnected vertices' triples in the networks, 
and the function C(k) restricted to classes of vertices with 
degree k allows to gather more detailed information. A 
possible weighted definition of the clustering coefficient 
is provided by the expression 



J (*) = 



1 



Si(ki - 1) 



E 



(Wij + w ih ) 



-Ijh, 



(9) 



where ctjh is equal to I if j and h are linked, and other- 
wise. The quantity c w (i) takes into account the weight of 
the two participating edges of the vertex i for each triple 
formed in the neighborhood of the vertex i; it measures 
the relative weight of the triangles in the neighborhood 
of a vertex i with respect to the vertex' strength |14j . C w 
and C w (k) are defined as the average over all nodes and 
over nodes of degree k, respectively. It is worth remark- 
ing that alternative definitions of the weighting scheme 
for clustering have been proposed in the literature |39|. 

Figure 0] displays for the North- American airport net- 
work the behavior of these various quantities as a func- 
tion of the degree. An essentially flat k nn (k) is obtained 
and a slight disassortative trend is observed at large k, 
due to the fact that large airports have in fact many 
intercontinental connections to other hubs which are lo- 
cated outside of North America and are not considered 
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FIG. 4: Assortativity and clustering for the North- American 
network. Circles correspond to topological quantities while 
squares are for affinity and weighted clustering. 



in this "regional" network. The clustering is very large 
and slightly decreasing at large k. This behavior is of- 
ten observed in complex networks and is here a direct 
consequence of the role of large airports that provide 
non-stop connections to different regions which are not 
interconnected among them. Moreover, weighted correla- 
tions are systematically larger than the topological ones, 
signaling that large weights are concentrated on links 
between large airports which form well inter-connected 
cliques (see also Qj| for more details). 



D. Betweenness Centrality 

A further characterization of the network is provided 
by considering quantities that takes into account the 
global topology of the network. For instance, the de- 
gree of a vertex is a local measure that gives a first in- 
dication of its centrality. However, a more global ap- 
proach is needed in order to characterize the real im- 
portance of various nodes. Indeed, some particular low- 
degree vertices may be essential because they provide 
connections between otherwise separated parts of the net- 
work. In order to take properly into account such ver- 
tices, the betweenness centrality (BC) is commonly used 
[40l lifl U2L l43j . The betweenness centrality of a node v 
is defined as 



E 



0~st 



(10) 



where o~ s t is the number of shortest paths going from s to 
t and a s t(v) is the number of shortest paths going from 
s to t and passing through v. This definition means that 
central nodes are part of more shortest paths within the 
network than peripheral nodes. Moreover, the between- 
ness centrality gives in transport networks an estimate 
of the traffic handled by the vertices, assuming that the 
number of shortest paths is a zero-th order approxima- 
tion to the frequency of use of a given node. It is generally 
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III. THE MODEL 
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FIG. 5: Scatter-plot of the betweenness centrality versus de- 
gree for nodes of the North- American air-transportation net- 
work. The red squares correspond to the average BC versus 
degree. 



useful to represent the average betweenness centrality for 
vertices of the same degree 



g(k) 



l 



E 9(v). 

v /k v —k 



(ii) 



For most networks, g(k) is strongly correlated with the 
degree k. In general, the largest the degree and the 
largest the centrality. For scale-free networks it has been 
shown that the centrality scales with k as 



g(k) ~ fc" 



(12) 



where /j, depends on the network [111 l42l Ejjj . For some 
networks however, the BC fluctuations around the be- 
havior given by Eq. 1)12(1 can be very large and "anoma- 
lies" can occur, in the sense that the variation of the 
centrality versus degree is not a monotonous function. 
Guimera and Amaral psl ] have shown that this is indeed 
the case for the air-transportation network. This is a 
very relevant observation in that very central cities may 
have a relatively low degree and vice versa. In Fig. [5] 
we report the average behavior along with the scattered 
plot of the betweenness versus degree of all airports of 
the North American network. Also in this case we find 
very large fluctuations with a behavior similar to those 
observed in Ref. |2^|. Interestingly, Guimera and Ama- 
ral have put forward a network model embedded in real 
space that considers geopolitical constraints. This model 
appears to reproduce the betweenness centrality features 
observed in the real network pointing out the importance 
of space as a relevant ingredient in the structure of net- 
works. In the following we focus on the interplay between 
spatial embedding, topology and weights in a simple gen- 
eral model for weighted networks in order to provide a 
modeling framework considering these three aspects at 
once. 



Early modeling of weighted networks just considered 
weight and topology as uncorrelated quantities [35| . This 
is not the case in real world networks where a complex in- 
terplay between the evolution of weights and topological 
growth does exist. For instance, if a new airline con- 
nection between two airports is created, it will generally 
provoke a modification of the existing traffic of both air- 
ports. In general, it will increase the traffic activity de- 
pending on the specific nature of the network and on the 
local dynamics. This effect is introduced in the model- 
ing of growing networks by the mechanism of a strength 
preferential attachment together with a dynamical redis- 
tribution of weights. Following this strategy it is possible 
to produce weighted networks with broad distributions of 
weights, connectivities, and stren gths , and correlations 
between weights and topology p4l25j. 

Here we consider a weighted growing network whose 
nodes are embedded in a two-dimensional space. As in 
the weighted model of Ref. [24| , it is reasonable to think 
that a newly created node n will establish links towards 
pre-existing nodes with heavy traffic or strength (hubs). 
Costs are however associated with distances and there is 
a trade-off between the need to reach a hub in a few hops 
and the connection costs. The cost naturally increases 
with the distance implying that the probability of estab- 
lishing a connection between the new node n and a given 
vertex i decays as a function of the increasing Euclidean 
distance d n i. As in the case of topological preferential 
attachment (i.e. connecting probability proportional to 
the degree 44] ), this trade-off can be expressed in two 
different ways: the connecting probability can decrease 
either as a power-law of the distance [43. Erl l47j or as 
an exponential with a finite typical scale 22] as it seems 
more natural for networks such as transportation net- 
works (see Fig. |2l or technological networks |4^. All the 
effects described in the next paragraphs are obtained in 
the case of an exponential decay exp(— d n i/r c ) but are 
also present in the case of a power-law (the effect 
of a decreasing scale r c is qualitatively the same as the 
effect of an increasing exponent a) . Eventually, the cre- 
ation of new edges will introduce new traffic which will 
trigger perturbations in the network. This model there- 
fore consists of two combined mechanisms: 

1 . Growth. We start with an initial seed of N vertices 
randomly located (with uniform distribution) on a 
2-dimcnsional disk (of radius L) and connected by 
links with assigned weight wq. At each time step, a 
new vertex n is placed on the disk at a randomly as- 
signed position x n (still according to a uniform dis- 
tribution). This new site is connected to m previ- 
ously existing vertices, choosing preferentially near- 
est sites with the largest strength. More precisely, 
a node i is chosen according to the probability 
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where r c is a typical scale and d n i is the Euclidean 
distance between n and i. This rule of strength 
driven preferential attachment with spatial selec- 
tion, generalizes the preferential attachment mech- 
anism driven by the strength to spatial networks. 
Here, new vertices connect more likely to vertices 
which correspond to the best interplay between Eu- 
clidean distance and strength. 

2. Weights dynamics. The weight of each new edge 
(n, i) is fixed to a given value wo (this value sets 
a scale so we can take wq = 1). The creation of 
this edge will perturb the existing interactions and 
we consider local perturbations for which only the 
weights between i and its neighbors j G V(i) are 
modified 

Wij -mjij+5-^-. (14) 

After the weights have been updated, the growth process 
is iterated by introducing a new vertex, i.e. going back to 
step (1.) until the desired size of the network is reached. 

The previous rules have simple physical and realistic 
interpretations. Equation (|13|l corresponds to the fact 
that new sites try to connect to existing vertices with 
the largest strength, with the constraint that the con- 
nection cannot be too costly. This adaptation of the rule 
"busy get busier" introduced in 0, |2{| allows to take 
into account physical constraints. The weights' dynam- 
ics Eq. I|14H expresses the perturbation created by the 
addition of the new node and link. It yields a global in- 
crease of wq + 6 for the strength of i, which will therefore 
become even more attractive for future nodes. 

The value of 6 characterizes the susceptibility of the 
network. If 5 < wq, the new link does not have a large 
influence. The case 6 « u>o corresponds to situations 
for which the new created traffic (on the new link n — i) 
is transferred onto the already existing connections in a 
"conservative" way. Finally, S > wq is an extreme case in 
which a new edge generates a sort of multiplicative effect 
that is bursting the weight or traffic on neighbors. 

The model contains two relevant parameters: the ra- 
tio between the typical scale and the size of the system 
77 = r c /L, and the ability to redistribute weights, 6. De- 
pending on the value of 77 and 6 we obtain different net- 
works whose limiting cases are summarized in Figure |H] 
More precisely, we expect: 

• For ?] > 1, the effect of distance is negligible and 
we recover the properties of the weighted model 
of Ref. |2^|. In this case, one obtains power-law 
distributions for connectivities and strengths with 
exponent 7 = (4<5 + 3)/ (26 + 1), as well as for the 
weights (exponent a = 2 + 1/6) . The strength and 
degree are linearly related by s w (k) ~ (26 + l)k. 
The effect of the redistribution parameter 6 is to 
broaden the various probability distributions, and 
to increase the correlations between topology and 
weights. Moreover, no correlations are introduced 
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FIG. 6: Different limiting regimes of the model depending 
on the value of its parameters 5 and n. 

between the topology of the network and the under- 
lying two-dimensional space, so that the distance 
strength s d grows simply linearly with the degree. 

• When 77 decreases, additional constraints appear 
and have consequences that we will investigate nu- 
merically in the following. Unless otherwise speci- 
fied, the simulations correspond to the parameters 
m = 3 (i.e. an average degree (k) = 6), and 6 = 1.0. 
We consider networks of size up to N = 10, 000, and 
the results are averaged over up to 100 realizations. 
All the observed dependences in 77 are essentially 
the same for other investigated values of 6. 

IV. NUMERICAL RESULTS 
A. Topology and weights 

At a purely topological level, the principal effect of a 
typical finite scale r c in the creation of new connections 
is to introduce a cut-off in the scale-free degree distri- 
butions In Fig. [7| we report the degree distribu- 
tion for a fixed value of 6 and decreasing values 77. A 
more pronounced cut-off appears at decreasing value of 
77 signalling the onset of a trade-off between the num- 
ber of connections and their cost in terms of Euclidean 
distance. The small-world properties of the network are 
as- well modified 22] : on the one hand, the increasing ten- 
dency to establish connections in the geographical neigh- 
borhood favors the formation of cliques and leads to an 
increase in the clustering coefficient (see inset of Fig.|SJ|. 
On the other hand, this same tendency leads to an in- 
crease in the diameter of the graph, measured as the aver- 
age shortest path distances between pairs of nodes. The 
diameter however still increases logarithmically with the 
size of the graph, as shown in Fig. 00 the constructed net- 
works do display the small- world property, even if strong 
geographical constraints are present. 

The correlations appearing between traffic and topol- 
ogy of the network are presented in Fig.[5]for two extreme 
cases of large and small 77. Strikingly, the effect of the 
spatial constraint is to increase both exponents (3 W and 
Pd to values larger than 1 and although the redistribu- 
tion of the weights [Eq. (|14fl ] is linear, non-linear rela- 
tions s w (k) and s d (k) as a function of k appear. For the 



7 




FIG. 7: Degree distribution P(k) for different values of 77 and 
5=1. The degree distribution is averaged over 50 networks 
of size N = 10 4 and minimum degree m = 3. 




FIG. 8: Average shortest path distances as a function of the 
network size N for different values of rj, with m = 3 and 5 = 1. 
Inset: clustering coefficient as a function of r\ for N = 10 4 , 
S — 1 (circles) and 8 — 5 (squares). The data are averaged 
over 50 networks. 



weight strength the effect is not very pronounced with 
an exponent of order (3 W w 1.1 for r\ = 0.01, while for the 
distance strength the non-linearity has an exponent of 
order (3d ~ 1-27 for rj = 0.02 (the value of the exponents 
(3 W and (3d depend on 77; see also 49] for a spatial model 
with (3d > 1). We show on Fig. the distance strength 
for two extreme situations for which spatial constraints 
are inexistent (77 = 10.0) or on the contrary very strong 
(r] = 0.02). 

The nonlinearity induced by the spatial structure can 
be explained by the following mechanism affecting the 
network growth. The increase of spatial constraints af- 
fects the trend to form global hubs, since long distance 
connections are less probable, and drives the topology 
towards the existence of "regional" hubs of smaller de- 
gree. The total traffic however is not changed with re- 
spect to the case 77 = 00, and is in fact directed towards 
these "regional" hubs. These medium-large degree ver- 
tices therefore carry a much larger traffic than they would 
do if global "hubs" were available, leading to a faster in- 
crease of the traffic as a function of the degree, eventually 




p d =1.27 



FIG. 9: Distance strength versus k for 77 = 0.02 and r\ = 10.0 
(the networks are obtained for 5 = 1, N = 10 4 , (k) = 6, and 
averaged over 100 configurations). When 77 is not small, space 
is irrelevant and there are no correlations between degree and 
space. When spatial effects are important (77 = 0.02 <C 1), 
non-linear correlations appear. We observe a crossover for 
k ~ 10 — 20 to a power-law behavior and the power-law fit 
over this range of values of k is shown (full lines). 



resulting in a super-linear behavior. Moreover, as previ- 
ously mentioned, the increase in distance costs implies 
that long range connections can be established only to- 
wards the hubs of the system: this effect naturally lead 
to a super-linear accumulation of s d (k) at larger degree 
values. 

Spatial constraints have also a strong effect on the 
correlations between neighboring nodes (Fig. I10[l . At 
large 77, a disassortative network is created, as is the case 
in most growing networks [5(J; as 77 decreases, k nn de- 
creases, and an increasing range of fiat k nn {k) appears: 
the tendency for small nodes to connect to hubs is con- 
trasted by the need to use small-range links. For small 
enough 77, a nearly neutral behavior more similar to what 
is actually observed in the airport network is reached. 
Moreover, the affinity of nodes to establish strong links 
to large nodes, measured by k™ n (k), goes from a flat be- 
havior at large 77 to a slightly assortative one at small 
77. In all cases, the weighted correlation fc™„(fc) remains 
clearly larger than the unweighted k nn (k), showing that 
links to busier nodes are typically stronger. 

A non-trivial clustering hierarchy is already displayed 
by the model without spatial constraints. As previously 
mentioned, the decrease of 77 leads to an increase of clus- 
tering. Moreover, the weighted clustering is always signif- 
icantly larger than the unweighted one, showing that the 
cliques carry typically an important traffic (see Fig. ITU)!. 
These effects are a general signature of spatial constraints 
as also observed in a non weighted network [2^ . 



B. Spatial constraints and betweenness centrality 

The spatial constraints act at both local and global 
level of the network structure by introducing a distance 
cost in the establishment of connections. It is therefore 
important to look at the effect of space in global topolog- 
ical quantities such as the betweenness centrality. The 
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FIG. 11: Betweenness centrality versus degree for differ- 
ent values of r\. Inset: Exponent \i (obtained by fitting the 
data for k > 10) of the betweenness centrality versus r\ (for 
N — 5,000, m — 3 and averaged over 50 configurations). For 
strong spatial constraints long-range shortcuts are very rare 
and hubs connect regions which are otherwise almost discon- 
nected which in turn implies a larger centrality of the hubs. 



FIG. 10: (a) Assortativity and (b) clustering versus k ob- 
tained for the model for 5 = 1, N = 10 4 , (k) = 6, 77 = 0.02 
(circles) and 77 = 1 (squares). Data are averaged over 100 con- 
figurations. Empty symbols refer to topological correlations 
while full symbols correspond to the weighted quantities k™ n 
andC"". 



betweenness centrality of a vertex is determined by its 
ability to provide a path between separated regions of 
the network. Hubs are natural crossroads for paths and 
it is natural to observe a marked correlation between g 
and k as expressed in the general relation g(k) ~ fc M . 
The exponent /j depends on the characteristics of the net- 
work and we expect this relation to be altered when spa- 
tial constraints become important. In the present model, 
Fig.llllclearlv shows that this correlation in fact increases 
when spatial constraints become large (i.e. when 77 de- 
creases). This can be understood simply by the fact that 
the probability to establish far-reaching short-cuts de- 
creases exponentially in Eq. I|13f) and only the large traffic 
of hubs can compensate this decay. Far-away geographi- 
cal regions can thus only be linked by edges connected to 
large degree vertices, which implies a more central role 
for these hubs. 

In order to better understand the effect of space on the 
properties of betweenness centrality, we have to explic- 
itly consider the geometry of the network along with the 
topology. In particular, we need to consider the role of 
the spatial position by introducing the spatial barycen- 
tcr of the network. Indeed, in the presence of a spa- 
tial structure, the centrality of nodes is correlated with 
their position with respect to the barycenter G, whose 
location is given by = ^ZjXj/iV. For a spatially or- 
dered network — the simplest case being a lattice embed- 
ded in a one-dimensional space — the shortest path be- 
tween two nodes is simply the Euclidean geodesic. In a 
limited region, for two points lying far away the probabil- 
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FIG. 12: (a) Betweenness centrality for the (one-dimensional) 
lattice case. The central nodes are close to the barycenter. (b) 
For a general graph, the central nodes are usually the ones 
with large degree. 



ity that the shortest path passes near the barycenter of all 
nodes is very large. In other words, this implies that the 
barycenter (and its neighbors) will have a large centrality. 
In a purely topological network with no underlying geog- 
raphy, this consideration does not apply anymore and the 
full randomness and the disordered small world structure 
are completely uncorrelated with the spatial position. It 
is worth remarking that the present argument applies in 
the absence of periodic boundary conditions that would 
destroy the geometrical ordering. This point is illustrated 
in Fig.Elin the simple case of a one-dimensional lattice. 

The present model defines an intermediate situation in 
that we have a random network with space constraints 
that introduces a local structure since short distance con- 
nections are favored. Shortcuts and long distance hops 
are present along with a spatial local structure that clus- 
ters spatially neighboring vertices. In Fig. ^| we plot the 
average distance d(G, C) between the barycenter G and 
the 10 most central nodes. As expected, as spatial con- 
straints become more important, the most central nodes 
get closer to the spatial barycenter of the network. 

Another effect observed when the spatial constraints 
become important are the large fluctuations of the BC. 



k/k 



FIG. 13: Average Euclidean distance between the barycenter 
G of all nodes and the 10 most central nodes (C) versus the 
parameter 77 (Here 5 = 0, N = 5, 000 and the results are 
averaged over 50 configurations). When space is important 
(ie. small n), the central nodes are closer to the gravity center. 
For large 77, space is irrelevant and the average distance tends 
to the value corresponding to a uniform distribution (r) un ij — 
2/3 (dotted line). 



FIG. 14: (a) Relative fluctuations of the betweenness central- 
ity versus k/k ma x for two values of 77 (N — 5,000 and the 
results are averaged over 50 configurations and binned). The 
fluctuations increase when 77 decreases (i.e. when spatial con- 
straints increase), (b) Number of anomalies N a (k) rescaled 
by the number of nodes Nj, versus k/k m ax for different values 
of n. The relative number of anomalies is larger when spatial 
constraints are large, especially for large k. 



Fig. 114b displays the relative fluctuation 



6g(k) = 



VWWj 

(g(k)) 



(15) 



where (Sg 2 (k)) is the variance of the BC and {g(k)} its 
average (computed for each value of k). The value of 
77 modifies the degree cut-off and in order to be able to 
compare the results for different values of r\ we rescale 
the abscissa by its maximum value k max . This plot 
(Fig. 114b ) clearly shows that the BC relative fluctuations 
increase as 77 decreases and become quite large. This 
means that nodes with small degree may have a rela- 
tively large BC (or the opposite), as observed in the air- 
transportation network (see Fig. [5] and H^|). In order to 
quantify these "anomalies" we compute the fluctuations 
of the betweenness centrality Afljv(fc) for a randomized 
network with the same degree distribution than the orig- 
inal network and constructed with the Molloy-Reed algo- 
rithm |5l| . We consider a node i as being "anomalous" 
if its betweenness centrality g(i) lies outside the interval 
[(g(k)) — aA^jv(A), (g{k)) + aAnjy(k)], where we choose 
a ~ 1.952 so that the considered interval would repre- 
sent 95% of the nodes in the case of Gaussian distributed 
centralities around the average. In Fig. 114b. we show the 
relative number of anomalies versus k/kkmax for differ- 
ent values of 77. This plot shows that the relative num- 
ber of anomalies N a (k)/Nk increases when the degree 
increases and more interestingly strongly increases when 
77 decreases. Note that since for increasing k the number 
of nodes Nk is getting small, the results become more 
noisy. 

The results of Figs. (|11I14|) can be summarized as fol- 
lows. In a purely topological growing network, central- 
ity is strongly correlated with degree since hubs have a 
natural ability to provide connections between otherwise 
separated regions or neighborhoods 01 . As spatial con- 
straints appear and become more important, two factors 
compete in determining the most central nodes: (i) on 



the one hand hubs become even more important in terms 
of centrality since only a large traffic can compensate for 
the cost of long-range connections which implies that the 
correlations between degree and centrality become thus 
even stronger; (ii) on the other hand, many paths go 
through the neighborhood of the barycenter, reinforcing 
the centrality of less-connected nodes that happen to be 
in the right place; this yields larger fluctuations of g and 
a larger number of "anomalies" . 

We finally note that these effects are not qualitatively 
affected by the weight structure and we observe the same 
behavior for 6 ~ or S ^ 0. 



V. CONCLUSIONS 

In this paper, we have presented a model of growing 
weighted networks introducing the effect of space and ge- 
ometry in the establishment of new connections. When 
spatial constraints appear, the effects on the network 
structure can be summarized as follows: 

• (i) Effect of spatial embedding on topology-traffic 
correlations 

Spatial constraints induce strong nonlinear correla- 
tions between topology and traffic. The reason for 
this behavior is that spatial constraints favor the 
formation of regional hubs and reinforces locally the 
preferential attachment, leading for a given degree 
to a larger strength than the one observed with- 
out spatial constraints. Moreover, long-distance 
links can connect only to hubs, which yields a value 
(3d > 1 for small enough 77. The existence of con- 
straints such as spatial distance selection induces 
some strong correlations between topology (degree) 
and non-topological quantities such as weights or 
distances. 

• (ii) Effect of space embedding on centrality 
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Spatial constraints also induce large betweenness 
centrality fluctuations. While hubs are usually very 
central, when space is important central nodes tend 
to get closer to the gravity center of all points. 
Correlations between spatial position and central- 
ity compete with the usual correlations between de- 
gree and centrality, leading to the observed large 
fluctuations of centrality at fixed degree. 

• (iii) Effect of space embedding on clustering and as- 
sortativity 

Spatial constraints implies that the tendency to 
connect to hubs is limited by the need to use small- 
range links. This explains the almost flat behav- 
ior observed for the assortativity. Connection costs 
also favor the formation of cliques between spatially 
close nodes and thus increase the clustering coeffi- 
cient. 

Including spatial effects in a simple model of weighted 
networks thus yields a large variety of behavior and in- 



teresting effects. This study sheds some light on the im- 
portance and effect of different ingredients such as spa- 
tial embedding or diversity of interaction weights in the 
structure of large complex networks and we believe that 
this attempt of a network typology could be useful in the 
understanding and modeling of real- world networks. 
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